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ABSTRACT 

We present an analysis of the deepest hard X-ray survey to date of about 2500 deg^ performed by the IBIS instrument on board 
INTEGRAL in the 20-60 keV band, with a total exposure time of 4 Ms. We find 34 candidate sources, for which we try to find 
counterparts at other wavelengths. The ratio of Seyfert 1 to Seyfert 2 is significantly more than the ratio found in the optical. This 
effect may be explained in the framework of the receding-torus model, but could also be due to absorption columns large enough to 
affect the 20-60 keV band. None of the predicted Compton-thick objects with lO^'* <Nn< 10^' cm"^ is detected unambiguously; when 
taking lower limits on A'h into account, the fraction of these objects is found to be lower than 24%. We do not see, but cannot exclude, 
a relationship between absorption and luminosity similar to what is seen in the 2-10 keV band. Our data suggests the possibility of a 
lack of objects with 10"' <Nii< 10^" cm"^, which could be expected if absorption has two origins, for instance a torus-like structure 
and the host galaxy. We find that the Log A'-Log 5 diagram of our sources is compatible with those obtained in other surveys in hard 
X-rays. Compared to models of the AGN population selected in the 2-10 keV band, the LogA'-Log5 diagram is generally in good 
agreement, but the A'h distribution is significantly different, with significantly less unabsorbed sources (A'h< 10^- cm"^) at a given 
flux limit compared to the models. In this survey, we resolve about 2.5% of the cosmic X-ray background in the 20-60 keV band. 
We also study the local hard X-ray luminosity function, which is compatible with what is found in other recent hard X-ray surveys. 
The characteristic luminosity Log ^'2o-60kcV ~ 43.66 is found to be a factor about 5 lower than the value observed in the 2-10 keV 
band. We find a space density of 10"^ AGN with L2o-60kev > 10*' per Mpc^^ and a corresponding luminosity density of 0.9 10^' erg 
s"' Mpc"\ 
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1. Introduction 

Recent progress in the understanding of the cosmological signif- 
icance of supermassive black holes stirred by both observations 
( [Magorrian et al. 1998 Ferrarese & Merritt 2000[l and numerical 



simulations (Kauffmann & Haehnelt"2000: Hopkins et al. 2005') 
prompts the making of a census of these objects as complete as 
possible. At distances too large to allow kinematics studies of 
stars in central regions of the galaxies, supermassive black holes 
are best revealed through their accretion of surrounding mate- 
rial. While these active galactic nuclei (AGN) are usually bright 
over a large fraction of the electromagnetic spectrum, the X-ray 
domain provides a privileged access to their population, if only 
for the reason that, at high galactic latitudes, the vast majority of 
X-ray sources are AGN. 

With the advent of powerful X-ray satellites, the 2-10 keV 
energy range is now very easily accessible. Both deep and wide 
surveys have been conducted and put together in order to study 
in detail the X-ray luminosity function of AGN up to cosmo- 
logical red shifts z~3 (|Ueda etal.|2003[|La Franca et al.|2005| l 
and higher ( Silverman et al. 2007| l, although sample sizes remain 
very small above z ~ 3. 

In spite of survey sensitivities in the /iCrab range, it is ex- 
pected that surveys conducted by XMM-Newton or Chandra 
may provide a significantly biased view of the AGN popula- 
tion against the most absorbed AGN. [Marshall et al.| ( |198 0) dis- 
covered with HEAO-1 an apparently diffuse X-ray emission at 
high galactic latitude. This so-called cosmic X-ray background 



presents a prominent peak around 30 keV, which has long been 
explained by the presence of highly absorbed AGN, with hy- 



^22, 



Compton-thick objects with A^h^ 10^"^ cm ^ ( 


Setti & Woltjer 


1989; Madau et al.||1994||Matt & Fabian||1994 


IComastri et al. 



1995). Even when they are intrinsically bright, these Compton- 
thick objects emit very little radiation below lOkeV and thus 
require deep X-ray observations. In a recent detailed modelling 
of the AGN population based on the known AGN population 
up to z ~ 3 and its extrapolation to higher redshifts, Gilli et aL] 
( |2007| l found that the population of Compton-thick AGN should 
be as large as that of moderately absorbed AGN. In fact a signif- 
icant fraction of local Seyfert 2 galaxies are found to be likely 
Compton-thick ( [Risaliti et al.|1999[|Guainazzi et al.|2005 1. 



Absorption is much less efficient in the hard X-ray domain 
(> 20keV) than in lower X-ray bands. Surveys in the hard X-ray 
have therefore the potential of detecting bright AGN with mini- 
mal bias in the A^h distiibution below ~ 10^^ cm"^. INTEGRAL 
and SWIFT are two satellites with such survey capabilities. 
Thanks to their large fields -of- view, the AGN population could 
be studied over the full sky with both INTEGRAL (Beckmann 
et al. 2006, S azonov et al |20 07) and SWIFT (Markwar dt eTaL 
2005; ,Tueller et al. 2007 1. Unfortunately, these surveys reach 
very limited sensitivities compared to their lower-energy coun- 
terparts. It is therefore very important to keep accumulating ex- 
posure time at high latitudes. In the case of INTEGRAL, which 
has a smaller field-of-view and a science program mostly geared 
towards the study of galactic sources, this means focusing on 



2 



S. Paltani et al.: A deep INTEGRAL hard X-ray survey of the 3C 273/Coma region 



a specific, relatively small high-latitude region. Provided a sus- 
tained observational effort is deployed over the next years of 
INTEGRAL operations these observations have however the po- 
tential of providing the deepest extragalactic hard X-ray survey 
for years to come, and is therefore an important goal. 

In this paper we study a ~ 2500 deg^ region of the sky cen- 
tered around 3C 273 and the Coma cluster with the goal of push- 
ing down as much as possible the sensitivity limit of hard X- 
ray surveys. This region is indeed the high-latitude region of 
the sky that has been the target of the deepest exposure with 
INTEGRAL. We study the population properties of the detected 
AGN, with a focus on their absorption properties. We compare 
the hard X-ray-selected AGN population with that resulting from 
medium (2-lOkeV) X-ray surveys. We also investigate the local 
AGN luminosity function in the hard X-rays. 

2. INTEGRAL observations of the 3C 273/Coma 
cluster region 

2.1. Data and processing 

Several INTEGRAL core-programme and open-time observa- 
tions have covered the sky region around Coma cluster and 
3C 273. We selected all available INTEGRAL pointings within 
30 degrees of a position located between these two sources, 
which resulted in 1660 pointings for a total elapsed observing 
time of 3 936 234 s and a dead-time corrected good exposure of 
2 733 202 s. Most pointings belong to four 5x5 dithering patterns 
repeated several times, plus a specific rectangular pattern used 
during the core-programme observation. 

Sky images in the 20-60 keV energy ranges were created 
from the data taken by the ISGRI detector of the IBIS imager 
on board INTEGRAL (Ubertini et al.|2003 ). Analysis has been 
performed with the Off-line Science Analysis software (OSA; 
[Courvoisier et al.||2003| p| version 6.0, using standard parame- 
ters. Good time intervals were built using a strong constraint on 
the attitude stability (deviation < 2 arcsec). Image cleaning was 
based on an input catalog of 34 sources with fixed source posi- 
tions and was applied independently of the source strength, al- 
lowing for negative source model to avoid introducing any bias 
in the process. 

To build the best possible mosaic of the field we excluded the 
outskirts of individual images, which are noisy and do not add 
much signal, as well as 10 individual images which had back- 
ground fluctuations larger than 1.1 in the significance map. 1650 
images were finally included in the mosaic, spanning revolutions 
from 0036 to 0464. 

The 3000x3000-pixel mosaic image was built in equatorial 
coordinates with a tangential projection using a factor-2 over- 
sampling when compared to the individual input sky images; this 
results in a pixel size of 2.4 arcmin in the center of the mosaic 
and of about 1.6 arcmin in the outskirts of the image, roughly 
40 degrees away from the center The photometric integrity and 
accurate astrometry are obtained by calculating the intersection 
between input and output pixels and weighting count rates ac- 
cording to the overlapping area. 

2.2. Properties of tiie mosaic 

Fig. [T] shows the exposure map of the field around the 
3C 273/Coma region. The total surface area is about 4900 deg^. 
The sky area having been exposed more than 10 ks is 2390 deg^, 

' http://isdc.unige.ch?Support+documents 




Fig. 1. Exposure map of the mosaic around 3C 273. The con- 
tours are located at 0, 10, 100, 300 and 500 ks respectively. 
Coordinates are right ascension and declination. Dots indicate 
the location of detected sources. 



while 1415 deg^ have been observed more than 100 ks. We point 
out that all quoted exposure times are effective exposures times, 
i.e. corrected for loss of efficiency due to dead time and off-axis 
observations. 

To investigate the quality of the mosaic, we study the 
distribution of the pixels' significance, which is expected to 
be Gaussian with mean and dispersion 1. Significant non- 
Gaussian tails are present when the full mosaic is studied. Such 
deviations are expected for sky areas having been observed a 
small number of times, where systematic effects have not been 
averaged out. This problem can be avoided by restricting the 
mosaic to the part where the effective exposure time is longer 
than lOks. Fig.[2]shows the significance distribution of the well- 
exposed pixels. We fitted a Gaussian distribution to the part of 
the histogram with cr < 3 to avoid the strong positive tail due to 
real sources. The negative branch of the significance histogram 
is very well fitted, with no evidence of excess in the tail. The 
sigma of the Gaussian is about 1.1, significantly larger than the 
expected value of 1 . The centroid of the Gaussian distribution is 
located at 0.009, which, while significant, has negligible impact 
on the significance levels of candidate sources. 

In view of the deviations found at low exposure times we 
shall completely discard the part of the mosaic with exposure 
shorter than lOks. The resulting mosaic contains almost 2 mil- 
lion pixels. 

3. Source extraction and selection 

3.1. Extraction method 

We build the candidate-source list from the mosaic by select- 
ing all pixels in the mosaic with a significance cTpixei larger than 
a given threshold small enough to ensure that every significant 
source has at least one pixel with a larger significance. We used 
o"pixei > 3. We use these pixels as starting points for the stan- 
dard OSA flux-extraction tool, mosaic_spec. This tool fits a 
Gaussian peak in a user-defined box centered around the input 
pixel. We used here a small 12xl2-pixel box to avoid that the 
fitted peak drifts towards nearby, more significant peaks; the box 
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Fig. 2. Histogram of the pixels' significance for the part of the 
mosaic exposed for more than lOks. The red curve is a Gaussian 
fit to the part of the histogram with cr < 3. 

size corresponds to 2.4 times the point-spread function's (PSF) 
FWHM in the center of the image (where pixels are the small- 
est). The Gaussian fit is performed letting the centroid free, but 
fixing the FWHM to that of PSF of ISGRI, i.e. 12 arcmin. As 
a result, we obtain for each starting-point pixel the coordinates 
of the peak with its associated uncertainty, a flux (again with its 
uncertainty), an exposure time and a revised significance cr, de- 
termined by mosaic_spec on the basis of the variance map. As 
the fits may converge several times on the same source, we used 
a distance-based algorithm to remove multiple detections of the 
same peak. 

3.2. Detection threshold 

Detection threshold obviously depends on the input mosaic, but 
it also depends on the exact method used to select and extract 
sources. In order to determine the detection threshold for the 
candidate sources, we investigate the distribution of the signifi- 
cance cr obtained if we apply the above extraction method start- 
ing from a random point on the mosaic. Fig.[3]shows the signif- 
icance distribution of the pseudo-sources obtained by drawing a 
large number of starting points at random over the part of the 
mosaic with more than lOks exposure. After removal of dupli- 
cates, we estimate that there are about 13 800 unique possible 
sources over the full mosaic. Their significance distribution fol- 
lows very well a Gaussian distribution, with a strong positive tail 
due to the fit occasionally converging on real sources. We fit the 
part of the distribution with significance cr < 4 with a Gaussian. 
The centroid of the Gaussian model is located at 1.38, and the 
variance is compatible with 1 . Thus, using the Gaussian approx- 
imation, one can derive the probability that source extraction at 
a fixed starting position returns a significance larger than a given 
significance ctq in absence of a source at this position. However, 
when looking for a priori unknown sources, we have to take into 
account that there are 13 800 unique sources in the mosaic. The 
probability to find a source with a significance cr > cry is there- 
fore given by the binomial distribution of the fixed-source prob- 
ability with 13 800 repetitions. The two probability cumulative 
distributions are shown in Fig. [3] The probability that a new can- 
didate source with a significance cr = 5 is real is about 20%. It in- 
creases to 85% for cr = 5.5 and to more than 98% for cr = 6. We 
point out that a number of candidate sources with significance 
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Fig. 3. Top: Significance distribution of the pseudo-sources ex- 
tracted randomly from the mosaic (blue line). The red line is a 
Gaussian fit to the part with cr < 4. The excesses in the sig- 
nificance distribution at cr > 4 indicate true sources. Bottom: 
Probability that the significance cr of a fake source extracted at 
a specific position is by chance larger than a given significance 
ctq (black line). The blue line gives the probability to find a fake 
source with significance cr > cro anywhere in the mosaic. 
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Fig. 4. Surface of the sky over which a given flux results in a 
significance cr > 5.5 (solid line) and cr > 5.0 (dashed line). The 
vertical dotted line shows the sensitivity limit of other hard X- 
ray surveys. 



just below the cr = 5 limit might turn out to be real sources. 
Using this probability distribution, we can now calculate the sky 
area over which a source with a given flux can be detected. An 
example is given in Fig.|4]for a significance threshold cr = 5.5. 
At the significance level cr = 5.5, the survey reaches a depth of 1 
count s"' over the totality of the 2390 deg^ with more than lOks 
exposure. 
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Table 1. Probability of chance coincidence for the candidate 
sources with counterparts in other catalogs. 



Object type 


Probability 


AGN(z<0.1) 


1% 


AON (all) 


3% 


X-ray sources 


5% 


Infrared sources 


12% 


Radio sources 


13% 


Galaxies 


26% 



3.3. Properties of the hard X-ray sources 
3.3.1. Source identification 

Tables |2] and [3] list the 34 candidate sources detected in the mo- 
saic with a significance cr > 5 together with their identifica- 
tions and basic properties, when available. Source identification 
is performed by searching for a counterpart in a radius depend- 
ing on the uncertainty on the source position. When the source 
cannot be obviously related to a known counterpart, we perform 
a tentative identification using both Simbad database|^from CDS 
(Centre de Donnees de Strasbourg) and the NASA Extragalactic 
Databas^ Counterparts in the error circles are selected based 
on the following criteria, in order of priority: (1) low-redshift 
AGN; (2) X-ray sources; (3) low-redshift galaxies; (4) Infrared 
or radio sources. Once a counterpart is chosen, class and red- 
shift information is taken from Simbad or NED, unless specified 
in Table |3] 

To quantify the probability of finding counterparts of a given 
type, we drew 100 random coordinates over the mosaic and 
searched for counterparts in a 3.5 arcmin radius using Simbad; 
this radius corresponds to the typical 90% position uncertainty 
of the least significant candidate sources. Table[T]gives the prob- 
abilities we found for different category sources. The probability 
of chance coincidence with an AGN is very low, especially if we 
restrict the redshift range to z < 0.1, which means that proba- 
bly all candidate sources identified with low-redshift AGN are 
true detections. We found however non-negligible probabilities 
of coincidence with apparently normal low-redshift galaxies, X- 
ray, infrared and radio counterparts, showing that such matches 
do not guarantee that the detection is not fake. In practice, we 
shall only consider that matches with low-redshift AGN indicate 
a secure detection. 



3.3.2. Flux and luminosity 

We calculate the source flux by fitting a standard AGN spec- 
trum to the 20-60 keV count rates extracted from the mosaic. We 
adopted a cut-off power-law with a slope T = 1 .9 and a cut-off 
energy = 200 keV (Pe rola et al.|[2U02| . We tested different 
choices of model parameters, for instance changing F by ±0. 1 
or setting Eq to lOOkeV or above 200 keV; the difference in flux 
has been found to be about 5%, which can be considered neg- 
ligible here. We find that a count rate of 1 s ' corresponds to a 
flux of 4.95 10"" ergs"' cm"^, which is about 5mCrab. When 
redshift information is available, luminosity is derived from the 
luminosity distance, using Hq - 70kms"' Mpc"', Qm - 0.3 and 
Qa = 0.7. 




Energy ( keV ) 

Fig. 5. Combined Chandra and ISGRI spectrum of 
MCG-03-34-064 in the 0.3-60 keV energy range. The 
source is well fitted by a model (solid line) consisting of a very 
absorbed power-law (dashed line) plus a reflection component 
(dot-dot-dashed line) and an Fe line (dotted line), a very steep 
component, probably due to scattered X-ray emission, is visible 
below 2keV (dot-dot-dot-dashed line). 



For comparison with other surveys in the hard X-rays, the 
20-60 keV flux of a source with a spectrum given by the 
above power-law is 1.55 higher than the 20-40 keV flux used 
in Beckmann et al.|(|2006); the 17-60 keV flux used in |Sazonov| 
et al. (2007 ) is 15% higher than the 20-60 keV flux; and the 14- 
195 keV (Tueller et al.,2007) is 2.16 times higher. 

Using these flux conversion factors, we find that previous 
hard X-ray surveys ( Beckmann et al.|2006 Sazonov et al.|2007 
Tueller et al. 12007) have all equivalent zero-area sensitivities in 



the 20 - 60keV domain of about 0.9 10"" erg s"' cm"^. At this 
flux limit, our survey still covers about 450 deg^, i.e. 15% of the 
full mosaic (see Fig. [4]i, while the ultimate depth of our mosaic 
is about 0.6 10"" erg s"' cm"^. On the other hand, these surveys 
cover a much larger fraction of the sky (> 75%) than ours (6%). 



^ http://simbad.u-strasbg.fr/simbad 

^ NED; ^http://nedwww.ipac.caltech.edu| 



3.3.3. Determination of the intrinsic absorption 

We searched the literature for previous X-ray measurements 
of the identified sources in the mosaic, in order to quantify 
their intrinsic hydrogen column densities A^h- For some objects, 
no absorption is detected beyond that originating from galac- 
tic hydrogen column density. In these cases, we list upper lim- 
its at the level of the galactic A^h- We note that, our field be- 
ing at high galactic latitude, these values are always on the 
order of 10^" cm"^. One of these objects, MCG-03-34-064, 
does not have a published spectrum, although it has been ob- 
served serendipitously by Chandra (Observation ID: 7373). We 
extracted and analyzed the spectrum, which turned out to be 
strongly absorbed, albeit still in the Compton-thin regime (A^h- 
0.74 lO^"* cm"^). The spectrum is presented in Fig.js] 

When no adequate X-ray measurement exists, we check 
whether there is a counterpart in the ROSAT all-sky survey 
bright source catalog (RASS-BSC; |Voges et alI][T9^ . This 
catalog is quasi-complete over the full sky down to a flux of 
0.05 ctss"' in the 0.1-2.4 keV band. Using the web version of 
HEASARC's PIMMS|^we calculated the expected count rate in 
the ROSAT band for different values of intrinsic A^h assuming a 



http : //heasarc . gsf c . nasa . gov/Tools/wSpimms . html 



S. Paltani et al.: A deep INTEGRAL hard X-ray survey of the 3C 273/Coma region 



5 



Table 2. Properties of the candidate sources detected in the mosaic. Coordinates are those of the excesses in the mosaic, e is the 
90% uncertainty on the position and 6 is the distance between the sources. Identification and classes are taken from Simbad or NED, 
unless specified in Table [3] cr is the source's significance. P is the probability that the source is real, not taking into account the 
identification. Star indicate a match with with low-redshift AGN or otherwise unquestionable AGN counterparts. 



Name 


RA 


Dec 




6 


Identification 


Class 


Exposure 


(T 


P 




J2000.0 


arcmin 






s 






IGR 112291+0203 


12 29 07 


+ 2 03 03 


0.20 


0.19 


3C 273 


Blazar 


699 100 


97.69 


1.00* 


IGR 112258+1240 


12 25 46 


+ 12 39 45 


0.35 


0.19 


NGC 4388 


Seyfert 2 


278 664 


79.90 


1.00* 


IGR 112106+3925 


12 10 34 


+39 24 39 


0.43 


0.38 


NGC4151 


Seyfert 1.5 


21 845 


42.77 


1.00* 


IGR 112396-0521 


12 39 37 


- 5 21 02 


1.64 


0.39 


NGC 4593 


Seyfert 1 


726431 


38.03 


1.00* 


IGR 112562-0547 


12 56 12 


- 5 46 33 


1.70 


0.87 


3C 279 


Blazar 


629 183 


10.98 


1.00* 


IGR J12392-1612 


12 39 10 


-16 11 41 


1.61 


1.29 


LED A 170194 


Seyfert 2 


249 172 


10.64 


1.00* 


IGR J12595+2755 


12 59 28 


+27 54 53 


2.07 


6.04 


Coma cluster 


Gal. cluster 


429 197 


10.42 


1.00 


IGR J 12233+0241 


12 23 21 


+ 2 41 05 


1.82 


0.92 


MRK50 


Seyfert 1 


636 864 


9.69 


1.00* 


IGR J12226+0414 


12 22 36 


+ 4 13 34 


1.96 


3.31 


4C 04.42 


Blazar 


590488 


9.68 


1.00* 


IGR J13092+1137 


13 09 10 


+11 37 26 


1.78 


1.31 


NGC 4992 


XBONG 


185410 


9.66 


1.00* 


IGR J12390-2719 


12 38 57 


-27 18 57 


2.06 


0.81 


MCG-04-30-007 


Seyfert 2 


34443 


8.98 


1.00* 


IGR J13383+0434 


13 38 17 


+ 4 33 51 


2.59 


1.32 


NGC 5252 


Seyfert 1.5 


58 992 


7.54 


1.00* 


IGR J12185+2948 


12 18 28 


+29 48 27 


2.52 


0.44 


NGC 4253 


Seyfert 1.5 


211 282 


6.71 


1.00* 


IGR J 13225 -1645 


13 22 29 


-16 44 41 


2.95 


1.51 


MCG-03-34-064 


Seyfert 1.8 


112 932 


6.55 


1.00* 


IGR J13041-0533 


13 04 05 


- 5 32 52 


2.77 


2.04 


NGC 4941 


Seyfert 2 


534 821 


6.13 


0.99* 


IGR J 12299+0305 


12 29 53 


+ 3 04 57 


3.67 


5.58 


IRXS J123013.6+030258 


X-ray source 


677 378 


5.92 


0.97 


IGR J 12069- 1448 


12 06 54 


-14 47 42 


3.59 


1.18 


2MASX J12065497- 1446335 


Galaxy 


147 802 


5.60 


0.89 


IGR J 13042- 1020 


13 04 13 


-10 19 58 


3.14 


0.60 


NGC 4939 


Seyfert 2 


414 987 


5.52 


0.84* 


IGR J12174-0131 


12 17 27 


- 1 31 19 


4.20 


2.64 


FIRST J121735.9-013001 


Radio source 


616797 


5.49 


0.83 


IGR J 120 11 +0649 


12 01 03 


+ 6 48 43 


3.23 


1.37 


LEDA 37894 


Seyfert 2 


286 255 


5.49 


0.82* 


IGR J13415+3023 


IJ 41 jZ 




J.JV 


A AQ 


MRK 268 


Seyfert 2 


189 232 


5.43 


0.78* 


IGR J12042-0756 


12 04 11 


- 7 55 44 


3.79 


1.62 


LEDA 157316 


Galaxy 


311 550 


5.38 


0.73 


IGR J12070+2535 


12 07 03 


+25 34 57 


3.42 


3.75 


IRAS 12046+2554 


Galaxy 


141 198 


5.37 


0.72 


IGR J13133-1109 


13 13 17 


-11 08 34 


3.59 


2.84 


IRXS J131305.9-110731 


Seyfert 1 


308 838 


5.36 


0.71* 


IGR J12060+3818 


12 06 02 


+38 17 48 


3.96 


2.80 


2MASX J12055104+3819308 


Infrared 


23 423 


5.26 


0.59 


IGR J13517-0042 


13 51 44 


- 42 27 


4.09 


3.40 


SDSS J135 133.23-004024.3 


Galaxy 


40371 


5.21 


0.52 


IGR J12130+0701 


12 13 01 


+ 7 01 22 


2.86 


1.03 


NGC 4180 


Seyfert 2 


420061 


5.20 


0.50* 


IGR J 11457- 1827 


11 45 41 


-18 27 29 


3.44 


0.26 


IH 1142-178 


Seyfert 1 


30405 


5.10 


0.36* 


IGR J12172+0710 


12 17 09 


+ 7 09 33 


3.31 


1.92 


NGC 4235 


Seyfert 1 


455 089 


5.09 


0.35* 


IGR J12136-0527 


12 13 37 


- 5 26 53 


3.84 


5.15 


IRXS J121353.4-053000 


X-ray source 


509473 


5.08 


0.33 


IGR J12310+1221 


12 31 00 


+12 21 28 


3.57 


3.26 


M87 


Radiogalaxy 


299 332 


5.08 


0.33* 


IGR Jl 1225-0419 


11 2231 


- 4 19 17 


5.07 


3.50 


2MASX Jl 1222455-0416096 


Galaxy 


21 144 


5.05 


0.28 


IGR J13353-1113 


13 35 17 


-11 13 10 


3.96 


6.00 


NVSS J133453-111300 


Radio source 


106440 


5.04 


0.27 


IGR Jl 1427+0854 


1 1 42 40 


+ 8 54 29 


3.92 


2.10 


2MASX J11424200+0852251 


Galaxy 


68 842 


5.03 


0.25 



power-law intrinsic emission with index F = 1 .9, taking into ac- 
count galactic absorption and source redshift. The method's ac- 
curacy is limited by effects of source variability, presence of soft- 
excess X-ray emission and extrapolation of a power-law whose 
actual index is unknown, but it is sufficient to obtain a moder- 
ately accurate estimate of A^h- Fig. |6] shows that RASS-BSC- 
derived A^h are estimated with an accuracy better than a factor 2. 
More importantly, it is able to identify correctly all the absorbed 
AGN(A?H> 10^2 cm-2). 

Absorption in the ROSAT 0. 1-2.4 keV band increases 
rapidly at hydrogen column densities around lO^^cm"^. 
Therefore the presence of a counterpart in RASS-BSC implies 
A^H^ 2 10^2 cm"2, depending quite moderately on the exact spec- 
tral shape. When no counterpart is present in RASS-BSC, we 
can constrain A^h by imposing that the expected count rate in 
the ROSAT energy band be below 0.05 cts s"', which translates 
into a lower limit A^h> 2 lO^^cm"^. For the objects for which 
we have no redshift, we have to make the further assumption 
that they lie at low redshift z « 1, otherwise this lower limit 
should be raised in their case. The case of 4C 04.42 shows that, 
should one of these sources have a redshift z ~ 1, the lower limit 
could be as high as lO^"' cm"^. Caution must however be taken 
for candidate sources without X-ray counterparts in the RASS- 



BSC, since fake detections have the highest chance of falling 
into this category. 

4. The population of hard X-ray selected AGN 

4.1. The mix of AGN types in hard X-ray selected samples 

In Table |2] we list all the sources in our observed area down to 
a significance level of cr = 5. As discussed earlier, because of 
the statistical properties of the mosaic, at this significance level 
there is a fairly large probability that the source is fake. As we are 
observing a field of high galactic latitude in a hard energy band, 
AGN are by far the most likely counterparts, with the exception 
of the Coma cluster ( Eckert et al.|2007 1. Conversely, association 
with a low-redshift AGN indicates that the candidate source is 
most probably real. From hereon in this section we consider only 
the 22 sources labeled with a star. 

Three blazars are very clearly detected, to which the possible 
identification of the radiogalaxy M 87 can be added. As 3C 273 
was the target of most of the observations, it is impossible to 
draw any statistics on the population of these objects. However, 
their redshift distribution is very different from that of the rest of 
the sources, as all blazars have redshifts higher than 0.1, reaching 
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Table 3. Properties of the candidate sources detected in the mosaic. Fluxes and luminosities are in the 20-60 keV energy range. 
niB is the apparent magnitude in filter B, unless specified in the footnote. A^h is the intrinsic hydrogen column density; values in 
parenthesis are derived from the RASS-BSC fluxes; lower limits mean non-detection in the ROS AT-BSC catalog. Unless mentioned 
in the comments, redshifts and class information have been obtained through Simbad or NED. Star indicate a match with with 
low-redshift AGN or otherwise unquestionable AGN counterparts. 



Name 


Flux 






Redshift 


login L„h^ 




10"" ergs 


cm-2 




10^2 cm-2 




erg s"' 


IGR J 1229 1+0203* 


11.56 ± 


0.12 


13.1 


< 0.01 


0.1583 


45.97 


IGR J12258+1240* 


15.27 ± 


0.19 


12.2 


26.92 


0.0084 


43.38 


IGR 112106+3925* 


30.74 ± 


0.72 


11.2 


7.5 


0.0033 


42.87 


IGR 112396-0521* 


4.48 ± 


0.12 


13.5 


< 0.01 


0.0090 


42.91 


IGR 112562-0547* 


1.37 ± 


0.12 


18.0 


0.05 


0.5362 


46.39 


IGR J12392- 1612* 


2.24 ± 


0.21 


10.9 


1.9 


0.036 


43.84 


IGR 112595+2755 


1.52 ± 


0.15 


13.5 


< 0.01 


0.0231 


43.28 


IGR 112233+0241* 


1.20 ± 


0.12 


15.5 


( 0.25) 


0.0234 


43.18 


IGR J 12226+04 14* 


1.23 ± 


0.13 


17.5 


( 7.1 ) 


0.9650 


47.08 


IGR 113092+1 137* 


2.36 ± 


0.24 


14.6 


90.0 


0.0251 


43.54 


IGR 112390-2719* 


5.48 ± 


0.61 


14.6 


10 


0.0250 


43.90 


IGR 113383+0434* 


3.25 ± 


0.43 


14.5 


2.89 


0.0230 


43.60 


IGR 112185+2948* 


1.45 ± 


0.22 


14 


< 0.02 


0.0129 


42.74 


IGR 113225-1645* 


2.02 ± 


0.31 


14 


73.85 


0.0165 


43.10 


IGR 113041-0533* 


0.84 ± 


0.14 


12 


44.98 


0.0037 


41.41 


IGR 112299+0305 


0.71 ± 


0.12 




( 2.5 ) 


- 


- 


IGR J 12069- 1448 


1.55 ± 


0.28 




> 2 


0.0184 


43.08 


IGR 1 13042- 1020* 


0.87 ± 


0.16 


1 1 


30.0 


0.0104 


42.33 


IGR 112174-0131 


0.70 ± 


0.13 


- 


> 2 


- 


- 


IGR 11201 1+0649* 


1.04 ± 


0.19 


15.3 


6.61 


0.0360 


43.51 


IGR J13415+3023* 


1 O^ _i_ 


r\ o^ 
V.Zi 


15.3 


> 2 


0.0399 


43.68 


Tr^D Ti on/10 ^7'^^^ 


1.00 ± 


0.19 




^ 






IGR 112070+2535 


1.42 ± 


0.27 


19.9 


> 2 


0.0477 


43.90 


IGR 113133-1 109* 


0.99 ± 


0.19 


15.6 


( 0.32) 


0.0344 


43.45 


IGR 112060+3818 


3.75 ± 


0.71 




> 2 






IGR 113517-0042 


2.70 ± 


0.52 


19.2 


> 2 


0.0531 


44.28 


IGR 112130+0701* 


0.81 ± 


0.16 


13.2 


> 2 


0.0070 


41.95 


IGRI11457-1827* 


3.29 ± 


0.65 


14.6 


(< 0.01 ) 


0.0329 


43.93 


IGR 112172+0710* 


0.76 ± 


0.15 


13.2 


0.15 


0.0080 


42.04 


IGR 112136-0527 


0.71 ± 


0.14 




( 2.0 ) 






IGR 112310+1221* 


0.94 ± 


0.19 


10.4 


0.02 


0.044 


43.65 


IGR 11 1225-0419 


3.54 ± 


0.70 


17.5 


> 2 


0.0535 


44.41 


IGRI13353-1113 


1.59 ± 


0.31 




> 2 






IGR 11 1427+0854 


2.08 ± 


0.41 


15.3' 


> 2 


0.0213 


43.34 


p band 



Comments 



1. 120O6 
()06| 



Shinozaki et al 



Nn:'Shuet al. ( 2007] 
A'H:Cappi etal. (2()66! 



Nh: 



Shinozaki et al. ( 2006^ 
Reeves & Turner 1 2000| 



Class: Masetti et al. ( 2()06c|; Nh: 
jVH: |Arnaud et al.^^2001j 



Sazonov et al. 



2005 



Class: Masetti et al. 



2006 



Class: Morelli et al. 
A'H:Shuet al. ( 2007 
Nh: Page etal. (2001 
A^h: This work (Chan dra): see text 



: Sazonov et 


il. 


po 


D5 


Tueller et al. 


2 


005 





Shinozaki et al. 



20061 



QSO 2MASS J12301552+0302546, z = 0.138, my = 17.5 ' 



A^H:'Bassani et al. 



1999 



No ob viously associated opt ical or infrared source 
Class: 



nij 



Masetti et al. 


2006b 


;Wh: 


Landi et al. 


2007 



14.1, niK = 13.0 



nij = 15.5, niK = 14.3 
Class: 



Masetti et al. 



2007 



Masetti et al. 



2006a 



Nh 

Sy 1 2MASX J12"1354 56-0530193. z = 0.066, my = 16.2 ? 
Nh^I 



Lieu et al. 



1996 



No obviously associated optical or infrared source 



almost z - 1, while the highest redshift for a Seyfert galaxy is 

0. 04. 

Out of the 18 remaining candidate sources, 17 could be iden- 
tified as Seyfert galaxies. Five of these Seyfert galaxies are clas- 
sified as Seyfert I's and eight as Seyfert 2's. The remaining 
Seyfert galaxies have intermediate types ( Osterbrock|| 1 98 1 [ l : 3 
Seyfert 1.5 and one Seyfert 1.8. The fraction of Seyfert galax- 
ies selected using optical spectroscopy has been determined by 
Hao et al.|p005) l in the Sloan Digital Sky Survey (SDSS). Out 
of 42435 galaxies, 1317 have an Ha full width at half maximum 
larger than 1200km s"', and should be classified as Seyfert I's, 

1. e. about 3%. Since the classification is based on Ha, this in- 
cludes intermediate types. Seyfert 2 galaxies are not so easily 
identified, and different selection criteria in the standard line- 
ratio diagrams ( [Veilleux & Osterbrock 1987[l lead to d ifferent 



fraction of Seyfert galaxies. Using the Kewley et al. ( 2001| i 



theoretical separation between star-forming galaxies and AGN, 
|Hao et al. (2005 ) find 3074 completely AGN-dominated narrow- 
line galaxies, i.e. about 7%. Using the more inclusive criteria 
of |Kauffmann et aT] ( |2003| l, which also selects AGN-starbursts 



composite galaxies, the figure increases up to 10700 galaxies, 
i.e. 25%. The ratio between Seyfert 1 and Seyfert 2 galaxies is 
therefore between 0.12 and 0.43. We find here a significantly 
larger ratio of 1. 13. Therefore a large fraction of Seyfert 2 galax- 
ies and/or Seyfert 2-starbursts composite are invisible in hard 
X-rays. It has been observed by Steffen et al. ( 2003 j l that the 
AGN 2-8 keV luminosity function was dominated by broad-line 
AGN at high X-ray luminosities and by narrow-line AGN at 
lower luminosities, contrarily to what is expected in the unified 
model of AGN ( | Antonucci 1 1 99 3 1 e.g.,). As we sample mostly 
the high-luminosity population, this can explain the discrepancy. 
IBarger e t al.| ()2005 ) analyzed in more detail the dependence of 
AGN type fraction on X-ray luminosity, and postulated that it 
could indicate that the covering fraction of the absorbing mate- 
rial decreases with luminosity, a scenario compatible with the 
receding-torus model used to account for a similar AGN-type 
dependence on radio luminosities ( Lawrence [1991 1. The limited 
statistics prevent us from checking whether this effect can com- 
pletely explain the fraction of broad-line AGN we observe here, 
but we point out that there could be an alternative explanation 
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20 21 22 

Log(N„) from RASS-BSC [ cm- 



24 



Fig. 6. Relationship between A^h measured in X-ray observa- 
tions and that derived from the counterpart in RASS-BSC. Coma 
and M87 are dominated by thermal emission in the soft X-ray 
band, and are therefore excluded. 



0.4 




0.1 



22 2A 
I,og(N„) [ ] 

Fig. 7. Histogram of the A^h distribution for the sources from 
Table [3] Sources are weighted by their probabilities. The blue 
histogram shows sources with measured intrinsic A^h's. The 
green histogram shows the sources whose A^h is derived from the 
presence of a counterpart in the RASS-BSC. The red histogram 
shows the sources without counterparts in the RASS-BSC, with 
logA^H evenly distributed between 22.25 and 25. 



to our result if a significant fraction of the Seyfert 2 galaxies are 
absorbed with column densities in excess of 10^^ cm"^, which 
would make them undetectable even in the hard X-rays. 

showed 



Deep X-ray surveys (e.g. Giacconi et al 



2001 



the existence of apparently normal galaxies with strong X-ray 
emission, the so-called X-ray bright optically normal galax- 
ies (XBONG; Fiore et aL| [2000). In this survey, we detect one 
XBONG, NGC 4992. While these sources consist most proba- 
bly of a mix of different kinds of objects ( jGeorgantopoulos ~&\ 
|Georgakakis| |2005 ), the high luminosity of this source (Lx ~ 
gj-g makes that it is quite probably a true AGN. 



4.2. Intrinsic absorption 

Figure |7] shows the distribution of the intrinsic hydrogen col- 
umn density A^h for the 34 sources detected in the mosaic and 



weighted by their probabilities. The fraction of absorbed ob- 
jects, i.e. those with A^h> 10^^ cm"^ is found to be 70%; if one 
discard the 11 sources without X-ray counterparts, this figure 
becomes 46%, making it a stringent lower limit. None of the 
23 sources with measured or estimated A^h are Compton-thick 



(A^H> 10 cm" ). We therefore obtain an upper limit to the frac- 
tion of Compton-thick objects of 24%. These figures are consis- 
tent with those found in previous hard X-ray surveys. In their 
modeling of the cosmic X-ray background, Gilli et al. ( 2007| l 
predict that, at the level of 10 " erg s ' cm"^, the fraction of ab- 
sorbed AGN (with A^H> 10^^ cm"^) should be 65% is perfectly 
compatible with our measurement. Their expected fraction of 
Compton-thick AGN (15%) is also compatible with our upper 
limit, although the lack of detection of any true Compton-thick 
object makes our constraint rather weak. 

There is a suggestion in Fig. |7] that there is a lack of ob- 
jects with A^H about 10^'"^^ cm"^. The chance of observing such 
a drop by chance from a flat distribution in Log A^h is about 
10%, which is not very significant, but the A^h distribution of 
SWIFT/BAT AGN (T ueller et aT] [2007)1 and in the previous 



INTEGRAL/IBIS survey of Beckmann et al. (2006) are both 



compatible with a similar drop. If verified on a larger sample, 
this drop could result from a double origin of the absorption: one 
absorber with large A'h, as in the torus model; and one absorber 
with moderate A^h in a galactic-type distribution. 



In the examination of local Seyfert 2 galaxies by Guainazzi 



et al. (20051, about 50% were found to be Compton-thick, i.e. 
with absorbing columns in excess of 10^"* cm"^. For the sample 
of 20-60 keV sources presented here, the median redshift (for 
the AGN with redshift information) is z = 0.023, and the me- 
dian 20-60 keV luminosity is lO'*^^ erg s"', i.e. we are probing a 
slightly more distant and higher luminosity population than that 
of the Guainazzi sample. 

We find that all Seyfert 2 galaxies are absorbed (A^h^ 
10^^ cm"^), while none of the Seyfert 1 galaxies are. Among the 
intermediate types, 2 are absorbed and 2 are not. This confirms 
the relationship between the absence of broad lines and obscura- 
tion found in local Seyfert galaxies (Cappi et al.'2006), which is 



again a natural consequence of the geometrical unification model 
of Seyfert galaxies (Antonucci 1993). 

The largest intrinsic absorption is actually observed for 
the XBONG NGC 4992, with A?h= 0.9 lO^"* cm'l This ob- 
ject seems similar to CXO U J031238.9 - 765134, which, ac- 
cording to Comastri et al.| |2002) , is a very strongly absorbed 
AGN. If Compton-thick, its hard X-ray luminosity should be 
about lO'^'^erg s \ close to that of NGC 4992.Using X-ray 
observations below lOkeV |Georgantopoulos & Georgakakis] 
(2005) found that XBONGs are a mix of quite different objects. 
However, many of these objects, like those presenting thermal- 
like emission with a temperature on the order of 1 keV, will not 
be detectable in the hard X-rays. Thus, it seems that hard X-ray 
selected XBONGs may to a large extent consist of a very ab- 
sorbed population of AGN. 

Figure |8] shows the relation between A^h and hard X-ray lu- 
minosity, as well as the fraction of absorbed sources in two lu- 
minosity bins, L < 10'*^ and L > 10*^ erg s"' respectively. A 
strong anticorrelation between the fraction of absorbed source 



and luminosity has been found in 2 - lOkeV surveys by Ueda 
et al. p003|) and later confirmed by several groups (|La Franca 



etal 



Dwe 



2005[|Georgakakis et al.|2006[|GiUi et al.|2007p , alfliougi 



ly & Page| ( |2006[ ) question this anticorrelation on the basis 



of extensive simulations, arguing that the discrepancy might re- 
sults from the combination cosmic variance (many results com- 
ing from the Chandra Deep Field South) and observation biases 



8 



S. Paltani et al.: A deep INTEGRAL hard X-ray survey of the 3C 273/Coma region 



20 - 



41 



A / 


\ / 




\' h 


l\ 




■ 










■ 

■ 






i 




■ 






1 ; 



0.8 



0.2 



42 



43 



45 



Log L^; 



[ erg 



4.3. Hard X-ray source counts 

Figure |9] top left shows the logA^-log5 diagram of the sources 
detected in the mosaic. To calculate this diagram we used all 
the candidate sources weighted by their probabilities of being 
real. 3C 273 and Coma cluster were however discarded from the 
log A^-log S diagram, because they were targeted by the obser- 
vations, and their presence in the mosaic is not due to chance; 
in addition. Coma cluster is a cluster of galaxies and its high- 
energy emission is most probably not due to the presence of an 
AGN ( [Eckert et al.|2007) . The faintest fluxes are approximately 
0.6 10"'' erg s~' cm~^. At this flux limit, the source surface den- 
sity is 0.033 + 0.010 deg-2. At the level of 10"" erg s"' cm^^, 
source surface density reaches 0.014 ± 0.0025 deg"^. The best-fit 
relationship is log = -1.55+0.13 log 5 - 1.86 + 0.07, 5 being 
the 20-60 keV flux in units of 10"" erg s"' cm"^. The slope is 
in quite good agreement with the expected -3/2 slope. 

We further point out that, at the sensitivity limit of this sur- 
vey, INTEGRAL resolves about 2.5% of the cosmic X-ray back- 
ground ( Gruber et al.|1999]|Churazov et al.|2007[l; this fraction 
is about twice that found in Beckman n et al.| ( |2006| l. A turnover 
in the log A^-log S relationship must be present at a flux between 



10-12 and lO-'-^erg s" 
medium X-rays. 



cm , similar to what is observed in the 



Be ckmann 
^ very 



4.3. L Comparison with otiier surveys 

At the level of 2 10"" erg s"' cm"^ between 20 and 60keV 
(which is a flux limit reached by all surveys) 

(2006 1 find a source density o f 4.3 10"^ deg 
close to our density of 4.7 10"^ deg"^. Tueller et al.| ( |2007| l, 
using SWIFT/BAT all-sky survey, find indeed a density of 
2.2 10"^ deg"2, lower than the INTEGRAL density by a factor 2. 
The authors point out that differences in Crab calibration might 
be the cause of this discrepancy, but the 15% difference in the 
Crab flux falls short of a factor 2 in source density. We checked 
that the assumption on the exact AGN spectral shape, for in- 
stance absorption, has a negligible eff'ect. Thus the discrepancy 
between INTEGRAL and SWIFT/BAT counts remains mostly 
unexplained. 

The 20-60 ke V source counts have a very similar slope and 
normalization to the bright end of the 2-10 keV source counts 



Fig. 8. A^H vs luminosity for all sources with significance cr > 5, 
the Coma cluster and blazars having been excluded. Blue filled 
squares correspond to objects with actual A^h measurements. 
Green open squares correspond to objects whose A^h's have been 
determined using their fluxes in RASS-BSC. Red arrows are 
lower limits based on the absence of X-ray counterparts. The 
heavy black line shows the fraction of absorbed AGN (A^h^ 
10^2 cm"2) in two luminosity bins, L < lO'*-' and L > lO'^^'erg 
s"'. The gray area shows the Icr uncertainties on these fractions. 



due to the limited bandwidth of Chandra compared to XMM- 
Newton. No correlation can be seen in our sample; the fraction of 
absorbed AGN is compatible with being constant around 70%. 
Among the three previous hard X-ray surveys, only |Sazonov] 
et al.| ( |2007] ) may have found an anticorrelation with some cred- 
ible significance. In any case, hard X-ray samples are still too 
small to exclude or confirm the anticorrelation claimed in the 
medium X-rays. 



From a compilation of 2-lOkeV surveys, Carrera et al. ( 2007[l 
find a best fitting slope of -1.58 and normalization of 485.3 deg" 



10 '"^ergs ' cm ^, 



equivalent to 0.01 13 deg at 



ig slope I 

at 1.17 X 

10"" erg s"' cm"^. Likewise, the 2-10 keV source counts pre- 
sented by |Moretti et al. ( 2003[ ) have a slope of -1.57 with 
a slightly lower normalization of 0.0082 deg"^ at S2-meV = 
10"") ergs"' cm~2. It should be pointed out that the 2-lOkeV 
source counts are only well defined at very faint fluxes (com- 
pared to hard X-ray surveys), and so a small uncertainty on the 
slope of the 2-lOkeV source counts can make a fairly large un- 
certainty when extrapolating up to ~ 10"" erg cm"' s"'. 

If we are seeing the same source population at 2-10 and at 
20-60 keV, then to reconcile the 2-lOkeV source count normal- 
izations of |Moretti et al.| ( |2003| ) and [Carrera et al. 1 (2007) with 
that seen here at 20-60 keV then the mean source flux ratio, 
'^ao-eokev/'^a-iokev. must be somewhere between 1 and 1.5. If 
we make the assumption that the average source spectrum is 
a power-law we require a mean slope of F ~ 1.5-1.8 to make 
the source count normalizations agree with each other Provided 
that sources without measured redshift do not lie at high red- 
shift (and hence they don't need any significant K-correction), 
the exponential cut-off at about 200 ke V ( Perola et al.|2002[ l in 
the average spectrum does not affect significantly the average 
spectral slope. This average slope of 1.5-1.8 is consistent with 
those found in AGN surveys at lower X-ray energies ( [PiconceUi 
et al.|2003]|Mateos et al.|2005 i. While this average slope, which 
is significantly harder than the typical F = 1 .9 intrinsic emission, 
implies the presence of absorbed sources, it suggests that at the 
20-60 keV flux level probed by our survey, there is no significant 
new population of Compton-thick sources that are not detected 
by surveys at lower energies, as the effective power-law slope 
between these two bands would be much harder. 



4.3.2. Comparison witli models of tine 2-1 keV AGN 
population 

In Fig. |9] top left we compare our source counts to the predic- 
tions of the 2-lOkeV population model of |Ueda et al.| ( [2003] l 
(hereafter U03). We predict the 20-60 keV source counts by in- 
tegrating the U03 model AGN population over the < z < 2, 
10"*' < i2-i0kev < 10'*^ erg s"' range. The conversion from rest- 
frame intrinsic (i.e. before absorption) 2-10 keV luminosity to 
observed frame 20-60 keV flux is made using the same spec- 
tral model adopted by U03, namely a power-law spectrum with 
F = 1.9, a cut-off rest-frame energy of 500 keV, and a reflec- 
tion component from cold material. With this spectral model, 
Sio-eokev = 1.24 X i2-iokev/47r<iL redshift ~ (^/l is the lu- 
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20-60 keV flux (erg s"' cm"^) 20-60 keV flux (erg s ' cm"^) 




Fig. 9. logA^-log5 diagrams for the sources detected by INTEGRAL. In all panels the solid lines are from our work and the 
dashed or dotted line are model extrapolations from U03. Top left: logA^-log5 for the full list of candidate sources with the Icr 
uncertainties. The solid line is the power-law best fit. Top right: logA^-log5 in two redshift bins; red: z < 0.05; blue: z > 0.05. 
Bottom left: logA^-logS' in two luminosity bins; red: Lio-eokev < lO"*"*'^; blue: L20-60keV > lO'*^'^. Bottom right: logA^-logS' in two 
A^H bins; red: Nh< 10-^ cm^^; blue: A^h> lO^^ cm^l 



minosity distance in cm). We use here the U03 model which 
includes a mix of unabsorbed and Compton-thin sources, but no 
Compton-thick sources; the Compton-thick population has in- 
deed not been measured and has been treated somewhat arbitrar- 
ily by adding a number of these objects equivalent to that of the 
Compton-thin ones. We have therefore assumed that absorption 
effects on this population in the 20-60 keV band are negligible. 

Under these assumptions, we can see that the U03 model 
provides a good match to both the slope and normalization of 
the total 20-60 keV source counts (see Fig. |9] top left), leaving 
little room for a significant additional population of moderately 
Compton-thick sources. We examine the population in more de- 
tail by splitting the sample into low (z < 0.05) and high red- 
shift sources, low (L2o-60keV < 10"*^^ erg s"') and high lumi- 
nosity sources, and low (A^h< 10^^ cm"^) and high absorption 
sources. In Fig.|9]top right we show the source counts and model 
predictions for redshifts either below, or above 0.05. Redshifts 
are unknown for 6 sources. On the basis of their optical mag- 
nitudes and alternative possible identifications we assume that 
all these sources, with the exception of IGR J12042-0756 and 
IGR J 12060+38 18, lie at z > 0.05 (and hence have 20-60 keV 
luminosities in excess of 10"^-^'^erg s"'). The observed counts 
and model predictions agree reasonably well given the rela- 



tively small numbers of observed sources. Above 52o-60keV ~ 
10"" erg cm"^ s"', the U03 model overpredicts the number of 
observed sources with z > 0.05. Irrespective of our redshift as- 
sumptions, the 20-60 keV integral source counts are dominated 
at all fluxes by low-redshift objects. 

In Fig. |9] bottom left we show the source counts and model 
predictions for observed sources above and below a luminos- 
ity of i20-60kev - lO^^*'^ erg s"'. For our given spectral model 
this corresponds to L2-iokeV = lO'*'*'*' erg s"' (as described 
above). The observed counts and model predictions are roughly 
in agreement, given the small number of sources. Note how- 
ever that the U03 model predicts roughly twice as many sources 
with high luminosity as with low luminosity, but that approxi- 
mately equal numbers of sources are observed above and below 
i20-60kev = lO'^^-^ergs"'. 

Absorbing column estimates or lower limits are available for 
our entire sample and so, in Fig. |9] bottom right, we show the 
source counts and model predictions separated into sources with 
A^H either greater than, or less than 10^^ cm^. Here it is clear that 
the U03 model is a poor predictor of the observed source counts. 
We see that the more absorbed sources constitute at least 2/3 of 
the total 20-60 ke V source counts over the flux range of the sam- 
ple, whereas the U03 model predicts equal number of absorbed 
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and non-absorbed sources over the luminosity and redshift range 
probed by the INTEGRAL observations. Taken at face value, it 
means that the observed and model Nh distributions differ signif- 
icantly. A number of our A^h measurements use indirect method 
based on the presence of the ROSAT flux or non-detection. A 
follow-up with more sensitive medium energy X-ray observa- 
tions is therefore needed to confirm this result. 

We have also investigated the dependence of the predicted 
source counts on the exact choice of spectral model. We find 
that the predicted source counts are relatively insensitive to the 
choice of cut-off energy, which is to be expected because the 
source counts are dominated by low redshift objects. The over- 
all normalization of the predicted source counts is rather depen- 
dent on the spectral slope and the size of the reflection com- 
ponent. A harder spectral slope or stronger reflection compo- 
nent increases the total 20-60 ke V counts predicted by the U03 
model. However, changing the spectral model has very little ef- 
fect on the predicted ratios of high to low redshift sources, high 
to low luminosity sources or high to low absorption sources. 

4.4. Hard X-ray luminosity function of local AGN 

We use the AGN firmly detected in this work (i.e. marked with a 
star in Table|2]to determine the AGN luminosity function (LF) in 
the 20-60 keV energy domain. However, because of the shallow 
depth of INTEGRAL surveys, the volume adequately sampled is 
quite smaU, making the high-luminosity tail of the LF, where ob- 
jects are very rare, impossible to determine. We therefore com- 
plement our sample with the sample of Beckmann et al. ( 2006 ) to 
study the part of the LF with L20-6O > lO** erg s"', using a clear 
cut in luminosity to avoid source duplication in these overlap- 
ping surveys. We convert the luminosities from the 20^0 keV 



to 20-60 keV domain using the flux conversion from Sect. 3.3.2 
We chose to study the redshift bin < z < 0.05 only to study the 
really local population and to minimize any effect of evolution; 
this has also the benefit of discarding the three blazars. We end 



up with 19 AGN from our sample and 2 from Beckmann et al. 
(|2006^. 

Fig. 10 left shows the luminosity function obtained from 
these samples. Because of the small number of sources, we use 
a parametric method to derive the LF. We assume the standard 
AGN luminosity function: 



3)' 



d<&(L) _ 
dlogL " (LIL*)-" + (LIL*)-P' 



(1) 



which describes a broken power-law, changing from index a to 
index yS at characteristic luminosity L* . From the area-detection 
threshold relationship, we build the probability function PdetiQ, 
which gives the probability that a source with a luminosity L and 
a redshift z < 0.05 is detected in this survey. Given a luminosity 
function <I>(L), we can derive the probability distribution of L 
in our survey. The parameters L*, a and /3 parameters are then 
determined using a maximum-likelihood (ML) test, based on the 
idea proposed by |Sandage et al.| ( |1979| l, modified to take into 
account our sensitivity map. As this method is independent of 
the normalization O*, we also estimate the LF using the standard 
non-parametric V/Vmax estimator O* is then obtained by fitting 
the y/Vmax LF estimates with the (S>{L) distribution from Eq. ([T]l, 
letting only O* free and fixing the other parameters to the ML 
values. 

To determine the uncertainties on the estimated parameters, 
we use the standard Gaussian approximation of the maximum- 
likelihood peak, i.e. the «-sigma confidence intervals for a are 



determined using the relationship: AX.(L* ,a,p) = exp(-Y), 
where X. is the likelihood and AX. is the ratio between the like- 
lihood for a given set of parameters and the maximum like- 
lihood. The (L*,a) uncertainty contours are shown in Fig. [TO] 
right. Confidence interval on <i>* is obtained by drawing sets of 
parameters {L*,a,/3) according to their likelihood and by fitting 
O* on the V/Vmax LF estimates. Finally, we obtain: 



logL* 
a 
P 



43.66 ^O:^^ erg s-i 

A OS +0.26 

a 19 +1.47 
•'■^^ -\.02 
1 1 9 +5.04 
^■^^-0.77 



(2) 



10"^ Mpc 



-3 



Fig. 



10 right shows, as an example, the uncertainty contours 
a plane. 3)* is badly constrained, because it is 



in the L* 

very strongly correlated to L*; therefore the LF normalization 
is known with a much better accuracy than the uncertainty on 
O* suggests. 



Figure 10 left also shows LFs determined from other sur- 
veys in sirnilar energy ranges ([Beckmann et al. 2006 Sazonov 
et al.|2007i[Tueller et al.|2007f | The LFs translated to the 20- 
60 keV range are very consistent with each other, except for the 
Bec kmann et al."(2006 ) LF, which shows a lower normalization 
by a factor about 2. The L* , a and fi parameters are consistent 
with each other in all studies. While the slopes a and fi obtained 
in the 2-lOkeV domain (e.g.. La Franca et al.|[2005 ') are quite 
comparable to ours, the characteristic luminosity L* is a factor 
about 4 higher than that observed here. As the conversion factor 
S 20-60 kev/'S 2-10 keV obtained by compai'ing the Log A^-Log S di- 
agrams was found to be between 1 to 1.5, this gives a difference 
by a factor about 5 in energy-corrected L* . This is formally com- 
patible with our study, but we note that all hard X-ray studies find 
L* values similar to ours, making the case for a statistical fluc- 
tuation rather weak. As |La Franca et alT] ( |2005[ ) use an evolution 
model over four redshift bins, the lowest one being < z < 0.5, 
it might be that the z = LF suffers from an extrapolation effect. 
It is however in line with the comparison of the source count di- 
agrams, where it appears that medium X-ray surveys see a less 
absorbed (and hence more luminous) AGN population. 

The spatial density of AGN with luminosity Lao-eokev > 
10"*' erg s"' is: 



0(L> 10"*') 



J41 



dO 
dlogL 



dlogL : 



1.03 lO^^Mpc"^ (3) 



This value is perfectly compatible with 'Sazonov eT al 
\200l\ (1.18 lO^^^Mpc-^) and with T ueller et al. ( 20071 
(1.23 lO^^Mpc"^), while |Beckmann et al. 12006 1 gives a lower 
value of 0.42 lO^^Mpc^l 



The luminosity density integrated over L 



20-60 keV 



WiL > 10^1) 



= /4>^ dlogL. 
= 0.90!|;;^^ lO^Vgs-i 



Mpc 



-3 



= 10"*' erg 



(4) 



Again , this is compatible with Sazonov et al. ( 2007| l and Tueller 



et al. 

T 



Mpc 



(|2007| (W(L > 10^1) = 1.02 10^^ and 1.03 10^^ erg 
Respectively), and about t wice the luminosity den - 



sity (0.41 10^** erg s ' Mpc ^ found in 
taking the flux correction into account. 



Beckmann et al. ( 2006 1, 



Tueller et al. 



j2007 I 



' We note that in the submitted version of the 
paper their LF is plotted with a wrong normalization, lower than the 
value they quote 
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Fig. 10. Left: AGN luminosity function in the 20-60 keV range. The points have been calculated using the V/Vmax estimator. The 
heavy solid line is the <1>(L) distribution from Eq. ([T]i with L* , a and j3 parameters determined using a maximum-likelihood method 
and obtained by fitting the V/Vmax points. Short-dashed line is the LF from Beckmann et al.| ( [2006 ); dotted line is that from 
Sazonov et al. ( 2007| l; long-dashed line is that from Tueller et al. ( 2007| l. Right: Uncertainty contours in the L*-a plane. The solid 
and dotted lines mark the 1,2,3 and 0.5,1.5,2.5 cr thresholds respectively. 



5. Summary and conclusion 

We presented an extragalactic survey of about 2500 deg^ 
with INTEGRAL/IBlS. With a flux limit of approximately 
0.6 10"" erg cm"^ s"' in the 20-60 keV range, this is the deepest 
hard X-ray extragalactic survey to date. We detected 34 candi- 
date sources, for which we searched for counterparts in other 
wavelength domains. We find that the ratio of detected Seyfert 1 
vs Seyfert 2 galaxies is larger than in optically selected samples. 
This suggests that current hard X-ray surveys are biased against 
Seyfert 2 galaxies, either because these surveys are limited to 
high-luminosity objects or because some sources are extremely 
absorbed sources, with A^h above 10^^ cm"^. 

We studied the distribution of absorption in our objects. 
While we do not find any Compton-thick object, using ob- 
ject without A^H determination and without counterpart in the 
ROSAT-BSC catalog, we find that at most 24% of the sources 
are Compton-thick. The fraction of absorbed object (Ah> 
10^^ cm"^) is between 46 and 70%. Although the small number 
of objects makes the result not significant, there's is a hint that 
the distribution could be bimodal, with absorption originating ei- 
ther from the host galaxy or from a thick structure with non-unity 
covering factor, like a dust torus. We do not find any relationship 
between A^h and luminosity, but, because of the small number 
of objects, we cannot exclude it. We note however that none of 
the current hard X-ray surveys find such relationship with any 
credible significance. 

The Log A-Log S diagram of our candidate sources has been 
found to reach 0.013 deg"' at the level of 10"" erg cm"^ s"'. We 
resolved approximately 2.5% of the cosmic X-ray background. 
Comparison with the 2-10 keV domain shows that the LogA- 
LogS' diagrams are compatible provided the average effective 
power-law index between the 2-10 keV and the 20-60 keV do- 
mains is r ~ 1.5 - 1.8, consistent with the 2-lOkeV slope found 



in medium X-ray surveys. This shows that there isn't any large 
population of bright Compton-thick objects missed in the 2- 
lOkeV surveys and appearing in the hard X-rays. 

We performed a detailed comparison of our source counts 
with those predicted from the population model of "Ued a et al.| 
(2003). Source counts are in general consistent, the biggest dif- 
ference being in the Log A-Log S diagram of sources separated 
in two Ah bins, Ah< 10^^ and Ah> 10^^ cm"^. As a large frac- 
tion of our Ah's were estimated based on ROSAT observations, 
our Ah distribution must be confirmed in follow-up observations 
in the medium X-rays. 

We present a truly local hard X-ray luminosity function 
(z < 0.05). We find a LF quite compatible with the latest all-sky 
INTEGRAL and SWIFT surveys. The characteristic luminosity 
L* is however a factor about 5 lower than that in the 2-lOkeV 
range. This discrepancy is seen in all hard X-ray studies of the 
local LF, and would imply, similarly to what is seen in source 
count diagrams, that the 2-10 keV population is less absorbed 
than that seen in the hard X-rays. 

The population of Compton-thick objects that is expected 
from models of the cosmic X-ray background and which should 
appear in the hard X-rays is still elusive. While the fraction of 
objects which are possibly Compton-thick is compatible with 
the models, the lack of clearly Compton-thick objects in our 
sample and the good match with the source counts from mod- 
els derived from medium X-ray observations make that the case 
for the existence of such population is rather weak. Follow-up 
observations of our sources without adequate Ah measurements 
are in progress, and may solve the puzzle. It seems however quite 
probable that we shall end-up with conflicting Ah distributions, 
which will have to be explained. 

While still suffering from low sensitivities compared to their 
modern counterparts working in the soft and medium X-rays, 
INTEGRAL and SWIFT are the only instruments currently 
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available to perform surveys above 15keV. The importance of 
this energy domain, which is unaffected by obscuration below 
~ 10^^ cm"^, is such that efforts to build statistically representa- 
tive samples must be pursued. 
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